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Abstract 



We present an output-sensitive algorithm for generating the whole set of fiats of a finite matroid. Given 
a procedure, P, that decides in Sp time steps if a set is independent, the time complexity of the algorithm 

f^ ' is 0{N^MSp), where TV and M are the input and output size, respectively. In the case of vectorial 

matroids, a specific algorithm is reported whose time complexity is equal to 0{N^Mcf), d being the rank 
of the matroid. In some cases this algorithm can provide an efficient method for computing zonotopes in 

^^ ' _ff-representation, given their representation in terms of Minkowski sum of known segments. 

^ ■ 1 Introduction 

A matroid is a structure, introduced by Whitney [T], providing an abstraction of the concept of independence 
that is common in different theories, such as linear algebra and graph theory. A finite matroid is defined 
as a pair (W,I), where W is a finite set, called ground set, and I is a collection of subsets of W, called 
independent sets, satisfying the following three properties [5]: 

^— s ' 1. The set I is not empty; 

^ ! 2. If A e Z and B C A, then Bel (hereditarity); 

C^ ' 3. If A, B G I and |A| < |B|, then there exists an element b £ B\A such that {b} U A G I {augmentation 

^— ^ , property), 

C| being the cardinality of set C. With a slight abuse of notation, hereafter we will denote by 6 U A the 

union of two sets A and {6}, the latter containing the single element b. In this paper we will refer often 

k>( \ to the concrete example of vectorial matroids, where the elements of W are vectors of a vector space and 

j_j ■ independent sets are the linearly independent subsets of W. 

Cd I We are interested in an efficient algorithm for computing fiats of a matroid. Flats are subset of W whose 

properties provide an alternative axiomatization of matroids. Ordered by inclusion, they form a geometric 
lattice [21 [3] ■ In the case of vectorial matroid, each fiat with maximal rank and properly contained in W 
(called hyperplane) can be associated with two facets of a zonotope up to translations. Zonotopes are polytopes, 
equivalently defined as Minkowski sum of segments or affine projections of cubes. They play an important 
role in several mathematical areas, such as hyperplane arrangements, box splines and partition functions [4]. 
They could turn to be useful also in some problems of quantum information. 

In general the number of flats grows exponentially in |W|. Thus, their computation has an exponential time 
complexity. However, in many practical problems the matroid has special properties that considerably reduce 
the output size with respect to the general case. For example, this occurs if the cardinality of some dependent 
subset of W is smaller than or equal to the matroid rank. An algorithm whose running time depends only 
on the input size does not take advantage of these special structures and its complexity is exponential in any 
case. Conversely, an output-sensitive algorithm, whose running time depends on both the input and output 
size, may require much less resources in cases of reduced output size. We will consider the extended notion 
of polynomial complexity that accounts for this output sensitivity. An algorithm is polynomial if its running 
time is polynomial in both the input and output size. 



In this paper we will present an output-sensitive algorithm for computing all the flats of a matroid. 
Assuming that there is a procedure that decides in Sp steps if a set is independent, the complexity of 
evaluating the flats is 0{N'^MSp), where M is the number of flats and N is the cardinality of W. The 
linearity in the number of flats is the significant feature that makes the algorithm output-sensitive. In the 
concrete case of vectorial matroid, the procedure P can be given for example by an algorithm that evaluates 
the rank of matrices. The overall complexity of evaluating the flats is 0{N'^M<P). We will also provide 
a specific optimization that requires an insignificant increase of computation space and reduces the time 
complexity to 0{N^M<P). Since the facets of a zonotope are identified by flats of a vectorial matroid up to 
translations, in some cases our algorithm can provide an efficient method for evaluating the _ff-rcpresentation 
of this particular polytope. The paper is organized as follows. In Sec. [5] we give some definitions, like rank, 
basis, closure and flat. We show that each flat can be represented through its bases. In sec. [31 we introduce a 
total order in the power set of the ground set W and use it to select a representative basis for each flat. Some 
properties of this basis are proved and then used in Sec. HI where we present the algorithm for generating the 
whole set of flats. In Sec. [SJ we show that our algorithm can provide in some cases an efficient method for 
evaluating the iJ-representation of a zonotope, given its representation in terms of Minkowski sum of known 
segments. 

2 Definitions and basic properties 

In order to define fiats, it is useful to introduce some concepts, such as basis, rank and closure. A basis is a 
maximal independent set, that is, an independent set that is not properly contained in an independent set. 
The last axiom of matroid theory implies that all the bases have the same cardinality, which is called rank of 
the matroid. Given a matroid (W,I) and any subset U of W, let J' be the collection of subsets of U that 
are in I, then the pair (U, J^) is a matroid. It is called the restriction of (W,I) to U '?|. Thus, each U C W 
has a rank, denoted by r(U), and a set of bases. 

The closure of a set U C W, indicated with cl(U), contains the elements w G W such that U and JJ Uw 
have the same rank, that is, 

cl(U) :=={we W|r(UUw) =r(U)}. (1) 

The closure operator is idempotent, that is, 

cl(cl(U)) =cl(U), VUCW. (2) 

A flat is a set F C W that is equal to its closure, that is, 

F is aflat ^F = cl(F). (3) 

A flat of rank n is called n-flat. The {d — l)-flat of a matroid of rank d is called hyperplane. Flats are 
analogous to vector subspaces. Indeed, in the case of vectorial matroids, it is possible to identify the flats with 
the subspaces linearly spanned by the vectors in the flats. In particular, n-flats correspond to ri-dimensional 
subspaces. 

By the idempotence property ^ we have that 

VX C W, cl(X) is a flat. (4) 

Trivially, the closure operator, with the sets of flats as codomain, is a surjective function, thus it is possible 
to represent every flat F through a subset whose closure is F. The minimal sets representing a flat F are 
independent sets with cardinality equal to r(F), as implied by the foUowings. 

Lemma 2.1 //X C Y and r(X) = r(Y), then cl(X) ^ cl(Y). 

Lemma 2.2 VX r(cl(X)) = r(X). 

As a direct consequence, we have the following lemma. 

Lemma 2.3 The sets with minimal cardinality that represent a flat F are all the bases o/F. 



Proof. Let S be a basis of F, thus, by definition of basis and rank, S C F and r(S) = r(F). By Lemma [2.11 
and definition of flat we have that cl(S) = cl(F) = F, that is, S represents the flat F. This set is also minimal 
because of Lemma [2121 Indeed, if R represents the flat F, then by definition cl(R) = F. By Lemma [2121 this 
implies that r(R) = r(F), that is, the cardinality of R is not smaller than the cardinality of a basis of F. 

Thus, flats can also be identified with the class of their bases, just as in linear algebra a set of k independent 
vectors identifies a fc-dimensional linear subspace. This provides a simplification in the representation of flats, 
since it is not necessary to enumerate the whole set of its elements. The mapping from the bases to the flats is 
surjective, but in general is not bijective. In the next section, we introduce a rule for selecting a representative 
basis that will turn to be fundamental for developing our algorithm. 

3 Denoting flats through representative bases 

To select a representative basis of a flat, we introduce a total order on the subsets of W and associate each 
flat with its first basis ("first" with respect to the total order). The order is defined as follows. 

Given a matroid (W,I), we order the elements of W by appending an integer i G [1,N] to each element 
Wi e W, N being the cardinality of W. Then we represent a subset U C W through a A^-digit binary number 
by setting the i-th digit equal to 1(0) if Wi is (not) an element of U, for every i e [1, N]. In other words, given 
any collection {6i, ..., 6fc} of indices with fe„-i < fe„, we label the subset X — {wb^ , wt^ , ■•■, wt^ } C W with the 
number 

k 

L{X) = ^2''--\ (5) 

n=l 

For example, 

L({u;i,W3,u;7}) = 10001012, 

which is equal to 69io in decimal basis. We call the most significant nonzero bit of a binary number b leading 
digit of b. 

By attaching the label L to each subset of W, we have introduced a total order in the power set of W. 
We associate each flat F with the basis that has the smallest label L, which we call pointer of F. In particular 
an n-pointer is the pointer of an n-flat. The pointer of F will be indicated with p{F). It is important to 
distinguish the pointer of a flat from its label L. The pointer of F is the label of its first basis, which in general 
is different from the label of F, unless F is independent. In binary representation, the number of non-zero 
digits of the pointer and label of a flat F is equal to r(F) and |F|, respectively. 

Theorem 3.1 Let pq be the pointer of a n-flat in binary representation, then the number s obtained from po 
by replacing the leading digit with is the pointer of an {n — l)-fiat. 

In order to prove it, we need the following property. 

Lemma 3.2 Let Ii and I2 be two independent subset 0/ W such that cl(Ii) — cl(l2), then, for every w G W, 
cl(Ii Uw) = cl(l2Uw). 

Proof. If ui e cl(Ii) — cl(l2), then the lemma is a direct consequence of the implication w e cl(X) ^ cl(X) = 
cl(X U w), which can be easily obtained from lemma [23] Thus, let us consider the case w ^ cl(Ii) = cl(l2), 
that is, we assume that IiUw and I2 U w are independent. We have to prove that if an element v ^ cl(Ii U w) , 
then V cl(l2 U w) and vice versa. Suppose that v ^ cl(Ii U w), then Ii U w U u is independent. Since also 
I2 U w is independent and |Ii U wUu| > II2 U w|, by axiom 3 of matroid theory (augmentation property) there 
is an element b ^ IiUwUv such that I2 U w U & is independent. Clearly, b is not in Ii, since cl(Ii) = cl(l2) 
and 6 ^ I2. Furthermore, b cannot be equal to w, thus b — v. This implies that I2 U it; U w is independent, 
that is, V ^ cl(l2 U w). Also the inverse implication is true. Thus, every element that is not in cl(Ii U w) is 
not in cl(l2 U w) and vice versa, that is, cl(Ii U w) and cl(l2 U w) are equal. D 

Proof of Theorem 13.11 Let Ii be the independent set with label s. Suppose that s is not a pointer, 
thus there is an independent set I2 with L{l2) < i(Ii) such that cl(l2) — cl(Ii). Let k be the position of 
the leading digit of po, then, by definition of s and Ii, Ii U Wk is the independent set pointed to by pq. By 
lemma 1321 both Ii U Wk and I2 U Wk are bases of the same fiat and furthermore i(l2 U Wk) < L{Ii U Wk), 



since i(l2) < ^(Ii), but this is impossible because pq = L{Ii U Wk) is a pointer. D 

This lemma implies that each (n+ l)-pointer can be generated from some n-pointer by setting one of the 
digits at the left of the leading digit equal to 1. For example, if a = IOOIOOII2 is a 4-pointer, then the number 
b = OOOIOOII2 is a 3-pointer for Theorem 13. II The 4-pointer a is generated from the 3-pointer b by replacing 
the 8-th zero digit of b with 1. Thus, given a collection of n-pointers, this procedure of replacement generates 
a set of labels that contains the set of all the {n + l)-pointers. In general the inclusion is strict, that is, the 
procedure of adding a bit 1 to an n-pointer does not necessarily give an {n + l)-pointer and we need a criterion 
for discarding labels that are not pointers. 

Theorem 3.3 Let 's ' and 'S ' be the label of an subset X and the position of the leading digit of s, respectively. 
Let Yi be the set obtained from X by removing the element Wj G X with j > i. The integer s is a pointer if 
and only if^ is independent and, for every wj. ^ X with k < S, w^. ^ cl(X) or w^. G cl(Yfc). 

Note that Y^ and Yj are not necessarily different ii i ^ j. In order to prove this lemma, we need three 
properties. The first one is known as the Mac Lane-Steinitz exchange property [2]. 

Lemma 3.4 Given a subsets C W and an element w G X, ifv G cl(X) andv ^ cl(X\u;), then c\(X\wUv) = 
cl(X). 

The second one is the basis exchange property [5]. 

Lemma 3.5 //Bi and B2 are two bases of a subset 0/ W and w G Bi, then there is an element v G B2\Bi 
such that Bi\w U v is a basis. 

Finally, the last property is a consequence of the augmentation axiom. 

Lemma 3.6 Let Y be a subset 0/ X G 1. If w ^ cl(Y), then there is an element v G X\Y such that 
w i cl(X\t;). 



Before proving theorem 13.31 let us prove the last property. 

Proof of lemma [3. 61 First, we assume that w G X. It is is clear that \iv = w then w ^ cl(X\w), since X 
is independent. Furthermore v G X\Y, since w G X and w ^ \Y, and the conclusion of the lemma is proved. 
Now we assume that u; ^ X. For the augmentation property, it is possible to construct an independent set 
by adding |X| — |Y U w] elements in X\Y to Y U w. The obtained set is equal to X U w minus some element 
V G X\Y. Thus, X\w U w is independent, that is, w ^ cl(X\w). D 

Proof of theorem 13.31 First we prove one direction of the implication and assume that s is the pointer 
of a fiat F. By definition X is independent. Suppose that the other part of the conclusion is false, thus there 
is an element Wfe ^ X with k < 6 such that Wk G cl(X) and Wk 4- cl(Yfe). Last condition and lemma [3T6l implv 
that there is an element wi G X\Yfc such that Wk ^ cl(X\w;). By definition of Y^ we have that I > k. Thus, 
since Wk G cl(X) and Wk ^ c\{X.\wi), by lemma |3^ we have that c\{X.\wi Uwk) — cl(X), that is, 'X.\wi U Wk 
is a basis of F, but this is impossible because L{'K\wi U Wk) < s {I is greater than k) and s is a pointer, thus 
one direction of the implication is proved. 

Let us prove the other direction. Suppose that s is not a pointer, then there is a number si < s such that 
i^^(si) = X and L~^{s) = X are bases of the same flat F. Let wi be the element with largest subscript I 
such that w; G X and wi ^ X. Since si < s, this element exists. By lemma [3751 there is an element Wk G X\X, 
such that ^\wi Uwk is independent, that is, Wk ^ c\(X\wi). The elements with subscript larger than I are in 
X if and only if are in X. Since Wk ^"^ and Wk G X, then k cannot be larger than I, thus we have that k < I. 
Because of this inequality, Yfc is a subset of X\w/. The relations Yfe C X\w; and Wk ^ c1(X\k;/) imply that 
Wk ^ cl(Yfe). It is also clear that Wk G cl(X), since Wk is in X and cl(X) = cl(X). D 

4 Algorithm for generating the flats 

Theorems 13. 11 and 13. 31 are the two key ingredients of our algorithm for calculating the fiats of a matroid. The 
idea is generating recursively the i-pointers from the lower-dimensional (j — l)-pointers. The flats are then 
generated from their pointers. More precisely, a set of labels is generated from each (z — l)-pointer by setting 



one of the digits at the left of the leading digit equal to 1. If Z is the position of the leading digit of an 
{i — l)-pointer and N is the number of elements in W, then a (i — l)-pointer generates TV — I labels. The set 
of labels generated from all the {i — l)-pointers contains the whole set of i-pointers. Theorem 13.31 provides an 
efhcient method for discarding labels that are not pointers. 

Since the structure of flats is unaffected by the presence of loops and parallel elements, we will assume 
without loss of generality that they are absent, that it, we assume that the matroid is simple j2]. Loops are 
dependent subsets of W with cardinality equal to 1. In the case of vectorial matroids, a loop is a zero vector. 
Parallel elements are pairwise dependent vectors. Denoting by N and d the cardinality and the rank of the 
matroid, respectively, the algorithm for generating the pointers is as follows. 

Input: The set of N 1-pointers {the elements in W} 

1: Set M equal to number of 1-pointers. {:— N} 

2: for i — 2,...,d - 1 do 

3: for j = 1,...,M do 

4; Set I equal to the position of the leading digit of the j-th {i — l)-pointer. 

5: ior S ^l + l,...,N do 

6: Generate from j-th (i — l)-pointer a label s by setting the S-th digit equal to 1. 

7: if s is a pointer {This is checked through theorem 13.31} then 

8: Store s as a new z-pointer. 

9: end if 

10: end for 

11: end for 

12: Set AI equal to the number of i-pointers. 
13: end for 
Output: The whole set of pointers 

The check at line[7]is performed through theorem l3.3l Thus, it requires to verify that the set X pointed to 
by s is independent and, for every w^ ^ X with k < S, the set XU w^ is independent or Y^ Uwk is dependent. 
Given a procedure, P, that decides with Sp steps if a set is independent, the time complexity of the check at 
line [7] is 0{NSp). The overall time complexity of the algorithm is 0{N'^MSp), where AI is the total number 
of flats, that is, M = J2iZi ^ij -^i being the number of i-flats. 

In the specific case of a vectorial matroid, the procedure P can be provided for example by a routine that 
evaluates the rank of matrices. Indeed any set X C W of a vectorial matroid can be seen as a d x |X| matrix 
(called column matroid), the columns being the elements of X. The rank of X is the rank of the matrix. 
The set is independent if the rank is equal to |X|. The time complexity of evaluating the flats by using this 
routine is O(N'^Md^). In this scheme the rank of the matrices is evaluated without taking advantage of the 
similarity of their structure. Indeed it is possible to reduce the time complexity by a slight increase of the 
computational space that exploits this similarity. Let po and Z be an (z — l)-pointer and its correponding 
independent set. / is the position of the leading digit of po- N ~ I labels are generated from po by setting the 
5-th digit equal to 1, where S is an integer ranging between I + 1 and N (see line [5] in the algorithm). We 
denote by X'^'^^ the sets associated with the generated labels. In line [7] of the algorithm, first we have to verify 
that X'^'^) is independent for each 6. Since sets with different values of S differ in one element and share the 
subset Z, the best strategy to decide if the sets X^^^ are independent is, first, reducing Z to the row echelon 
form using row operations and, then, performing the same raw operations on the added column in each X' ' . 
The corresponding time complexity is 0{d^N), taking into account that d < N. A similar strategy can be 
used also in checking the independence of X U Wfc and the dependence of Y^ U Wk- In this way it is possible 
to reduce the overall complexity of generating the flats to 0{N^Md^). 

5 Minkowski sum of segments: zonotope 

The algorithm for the computation of flats can be useful in some cases for calculating the _ff-representation 
of a zonotope when it is represented as Minkowski sum of known segments. The Minkowski sum of two sets 



A and i? in a vector space is the set obtained by adding every vector of A to every vector of B, that is, 

A + B = {d+b\ae A,b(EB}. (6) 

The zonotope is a polytope defined as the Minkowski sum of segments. Up to a translation, it is the set of 
vectors 



M 
V 

k=0 



= '^XkWk, (7) 



where A^ G [0:1]. Each vector Wk and the zero vector are the two vertices of each segment summed up. 
The set W = {wk\k € [1 : M]} is the ground set of a vectorial matroid, whose independent sets are the sets 
of linearly independent vectors. Let d be the dimension of the zonotope, that is, the maximal number of 
independent vectors in W. Without loss of generality, we assume that d is also the dimension of the vector 
space. 

Each facet is parallel to a (d — l)-flat of the matroid, thus its normal vector is orthogonal to any basis of 
the {d — l)-flat. The computation of each normal vector by a set of d — 1 vectors has a complexity that scales 
like d'^ . We denote a normal vector with Hi, where the subscript i is an integer that goes from 1 to M^-i, 
Md~i being the number of (d— l)-flats. It can be proved that for each vector Hi there exist two facets defined 
by the inequalities 

Hi ■ X <^^ 0{Hi ■ Wk)ni ■ Wk (8) 

k 

and 

Hi- x> -'^9[-Hi-Wk)ni-Wk, (9) 

k 

where 9{x) is the step function 6{x > 0) = 1, 9{x < 0) = 0. These inequalities define the zonotope in 
-ff-representation. Thus, the computation of the zonotope in iJ- representation is achieved by evaluating the 
(d — l)-flats of the matroid W. 

The output in this problem is the set of half-planes, thus its size is Md-i- In general our algorithm for the 
computation of the flats does not allows us to solve this problem in polynomial time with respect to the output 
size, since the algorithm presented in the previous section is linear in the total number of flats M and in the 
worse case M could be exponentially greater than Af^-i. However, in many practical problems M can scale 
linearly in Md~i and the input size. Suppose for example that the vectors Wk are in general position and the 
matroid rank d is smaller than N/2. The number of /c-flats is k\iN~k)\ ^^'^ grows monotonically in k{< d—1). 
This implies that M scales at most like dM^-i- This linear scaling can be present also in the case of special 
structure for which an output-sensitive algorithm provides an advantage. It is worthwhile to note that the 
best algorithm for the evaluation of a zonotope in ^-representation has a complexity that is quadratic in the 
output size in any case [5] . In the subclass of problems where the number of overall flats is a linear function 
of the number of hyperplanes, our method for the generation of the zonotope in ^-representation is linear in 
the output size. An open question is determining how much large is this subclass 
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